Non-exponential relaxations in disordered conductors 
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We show that, in low dimensional conductors, the quasiparticle decay and the relaxation of the 
phase are not exponential processes. In quasi-one dimension, they scale as e { ' N > where 
the characteristic time n n , identical for both processes, is a power T 2 ^ 3 of the temperature. 
This result implies a distribution of relaxation times. 



The problem of dephasing in the presence of electron-electron interaction in disordered con- 
ductors was first addressed in the pioneering work of Altshuler and Aronov 1 . It has been revived 
by a series of recent experiments 2,3,4 . Another very related question is the understanding of the 
time evolution of a quasiparticle state which governs the relaxation towards equilibrium. This 
quasiparticle decay is usually described by a characteristic time, implying implicitly an expo- 
nential process. Here we show that in low dimension, and in particular in quasi-one-dimensional 
disordered wires, this decay is faster than exponential. We focus on two aspects, the decay of 
a quasiparticle state and the loss of phase coherence. It has been shown that both processes 
are related 5 . They are described by the same characteristic time scale. Here we show that both 
processes are non exponential 6, . 

We start by considering the lifetime of a quasiparticle state, initially in a state of energy e. 
Its interaction with another quasiparticle of energy e' leads to final states at energies e — ui and 
e' + <jj. At zero temperature, the lifetime is given by 
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where uq is the density of states per spin direction. The kernel K(iv) is the disordered averaged 
squared matrix element of the screened Coulomb interaction. It depends only on the energy 
transfer ui and it can be written as 1 
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In this expression P(q,tu) = l/(—iuj + Dq 2 ) is the classical diffusion probability, U(q,uj) is 
the dynamically screened coulomb potential U(q,ui) = (l/2po)Dq 2 /(— iuj + Dq 2 ) and po is the 
density of states per spin direction and per volume unit. The kernel can also be written under 
the form 
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so that the lifetime reduces to 
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Replacing the discrete sum over q by an integral leads to 



where ct\ = V2/ir 2 , a.2 = l/2vr 2 , a 3 = V2/2ir 3 and E c = D/L 2 is the Thouless energy. From 
Eq. (1), one obtains the well-known 1 power-law dependence of the lifetime l/r ee (e) oc e d l 2 . 

At first sight, one could consider that, at finite temperature, the quasiparticle relaxation 
rate is given by l/r ee (e, T) oc max(e <i / 2 , T d l 2 ). This is not correct for d < 2 because of the 
importance of the processes with low energy transfer 1 . However, it has been found that the 
quasiparticle relaxation is still exponential with a characteristic time which varies like T 2 / 3 in 
one dimension 1 . Here we show that in low dimensions d < 2, the relaxation is not exponential. 
To that purpose, we start from the expression 
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that generalizes Eq. (1) to finite temperatures. The thermal function F(e,e',u) depends on the 
Fermi factors f e = l/(e l3e + 1) and has the form 

F(e,e',U J ) = f e ,(l-fe-u J )(l-fe' + , J ) + (l-fe>)fe-u J fe^ • (7) 

The first term is larger for e > 0, it describes the relaxation of an electronic state above the 
Fermi level. The second term, larger for e < 0, accounts for the relaxation of a hole state 
within the Fermi sea. The two terms are equal for e = 0. Integrating over e', and after some 
straightforward manipulations, it rewrites 
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This expression can also be obtained from the imaginary part of the self-energy, a calculation that 
incorporates implicitly both processes involved in Eq. (7) 8 . r ee (e,T) is defined as the lifetime 
of a quasiparticle. It can also be interpreted as being the characteristic time that describes the 
relaxation towards equilibrium. Indeed, it can be defined from a Boltzmann equation for the 
energy distribution n e 1,9 
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The relaxation term contains two contributions that respectively describe the quasiparticles 
leaving a given quantum state ( "out" term) and incoming in this state ("in"). At equilibrium 
n e equals the Fermi factor f e = l/(e /3e + 1) and the relaxation term is zero. Linearising around 
the equilibrium distribution n e = f e + 5n t , one obtains a relaxation of the form dSn e /dt = 
—Sn e /T ee (e,T) with the same characteristic time given by (6). The lifetime of a quasiparticle 
state can also be interpreted as the relaxation time of the energy distribution. 

We now consider the lifetime of a quasiparticle at the Fermi level and we denote it by Tj„(T). 
From Eq. (8), it is given by 10 
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In d = 3, the integral converges, leading to the well-known T 3 / 2 behavior 1 . For d < 2, the 
integral diverges at low energy transfer. It has been argued that the integral has to be cut-off 
at low energy: since the lifetime of a quasiparticle is finite, no energy transfer can be smaller 
than this inverse lifetime 1 . Consequently, the lifetime is solution of the self-consistent equation 
(for simplicity the thermal factor in Eq. (10) has been replaced by an upper cut-off at to ~ T) : 



(12) 



This argument gives a characteristic time Tj n (T) which scales as T 2 / 3 in one dimension : 
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A = 1/vo is the interlevel energy spacing. 

We argue that the divergence at low energy transfer in (11) has a deeper significance 6 . It 
is the signature of a non- exponential relaxation. The reason goes as follows: the quasiparticle 
state decay cannot be exponential since, after a time t, the energy transfer cannot be defined 
with a precision better than \ jt (Heisenberg inequality) . Thus the decay rate cannot be constant 
in time 6 . Instead of assuming a relaxation of the inital state of the form V(t) = e~ l l Tin with eq. 
(11), we find that the probability V(t) for the quasiparticle to stay in its initial state is given by 
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an expression valid for times t S> 1/T. In one dimension, one has 7 : 
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which leads to a non-exponential behavior for the relaxation of a quasiparticle state. 

V(t,T) ~ e -['M»(T)]3/ 2 d=1 (15) 

with the same characteristic time as in Eq. (12). Similarly, in two dimensions, starting from 
(13), one finds a logarithmic correction to the exponential decrease 



V(t, T) ~ e T in ln « d = 2 (16) 



where T in {T) oc AT/E C . 



Now, one can wonder whether this peculiar behavior of the relaxation of a quasiparticle 
state appears in the time dependence of other quantities. To that purpose, we consider the 
phase coherence time. It is defined as the lifetime of the Cooperon in the presence of other 
electrons. The Cooperon is the quantum correction to the return probability that consists in 
pairs of time reversed trajectories. Altshuler, Aronov and Khmelnitskii 12 have shown that the 
effect of other electrons can be accounted for by a fluctuating electric potential V(r,r) whose 
characteristics are given by the fluctuation-dissipation theorem: 

The Cooperon contribution to the return probability takes the form 

P c (r,r,t)=P^(r,r,t)^ r ^) Tc (18) 

where is the cooperon in the absence of the fluctuating potential. = Q(r, t) is the relative 
phase for a pair of time reversed trajectories at time t : 
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One of the trajectories propagates from time r = to r = t, while the time reversed trajectory 
propagates from r = t to r = 0. We define f = t — r and (• • is the average taken both 
on the distribution of diffusion trajectories ((• • -}c) and on the thermal fluctuations ((• • -)t) of 
the electric potential. The thermal fluctuations are gaussian so that the thermal average (e**) T 
satisfies : 

( e ^) T = e -H* 2 ) T . (20) 
Let us start with the calculation of (& 2 ) T - Inserting (17) into (19), one finds : 
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In one dimension, the phase fluctuation becomes 

= [w-rmW . (22) 

It depends on the trajectory r(r) and we still need to calculate the average 

(e-^ 2 )^ (23) 
over the distribution of diffusive trajectories. A first approximation consists in assuming that 



e 2 \ itj = e 2 ^ ' T > c . Averaging over the diffusive trajectories, we obtain that |r(r) — r(r) 
scales as \fl so that, from (22), one has 
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where we have defined the characteristic time 
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so that we find, within this approximation : 
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A similar behavior has been found in 11 . Apart from a numerical factor, the characteristic time 
(25) is the same as the inelastic time (12) for the decay of a quasiparticle state. Moreover Eq. 
(24) shows that the phase relaxation has the same temporal dependence as the quasiparticle 
relaxation (15). 

However, the result (26) is not yet fully correct since we have made an approximation in 
replacing the average (• • •) of the exponential by the exponential of the average. The exact 
behavior of the phase relaxation (23) can be obtained from a functional integral approach 12 
that shows that the Laplace transform 

P 7 (r,r) = J dtP c (r,r,t)e-^ 

is given by : 
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with Ai et Ai' being respectively the Airy function and its derivative 13 and r 7 = I/7. The 
probability P c (r, r, t) could be obtained from inverse Laplace transform of Eq. (27). Its time de- 
pendence is clearly non exponential. Here we obtain the time dependence of the phase relaxation. 
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Figure 1: Behaviour of (e l ^ t ')T,c- The continuous line is the exact result (29). The dotted line is obtained from 
the small time expansion (26). The dashed line shows the exponential fit e ~ t//2Ti ™ . 

It is of the form (e**) TC = /(£/ t at)- Since in one dimension, one has P c (0) (r,r,t) = l/y/torDi, 
the unknown function /(t/Tjv) is solution of 

Jo V^Di \tnJ 2 V D Ai' (Tjv/r 7 ) v ; 

The inverse Laplace transform is obtained by noticing that both the Airy function and its 
derivative are analytic and non meromorphic functions whose zeroes lie on the negative real axis. 
Then, by performing the integral in the complex plane with the residues Res(e s *Ai(s)/Ai'(s)) = 
e~'"™'*/l u nl where the u n are the zeros of Ai'(s) given at a very good approximation by \u n \ = 
(^(n — |)) 2//3 13 , we obtain the analytic function 

(e l %C = J^Y,- 1 ^T- (29) 



At small times t < Ti n , it behaves like eq. (26). At large time, the relaxation is driven by 
the first zero of the Ai' function, namely (e**) TC — y^nt /Ti n e~\ Ul \ t / Tin /\u±\ with \m\ ~ 1.019. 
Clearly, the relaxation (29) is never exponential. It appears as a distribution of relaxation times 
Tin/\un\ which is at the origin of the rather unexpected compressed exponential behavior of the 
quasiparticle decay and of the Cooperon phase relaxation. 

In conclusion, we have shown that the decay of a quasiparticle state in a low dimensional 
disordered conductor is not exponential. In quasi-ld, it is of the form e~(*/ Tin ) 3/2 . We have also 
calculated the relaxation of the phase in a quasi-ld conductor and we have shown that is it 
also described by the same time dependence as for the relaxation of the energy, with the same 
characteristic time. 
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